The instability of Alexander-McTague crystals and its implication 

for nucleation 



W. Klein 

Center for Non-Linear Studies, Los Alamos National Laboratory, Los Alamos, NM 87545 
Department of Physics and Center for Computational Science^ 
Boston University, Boston, MA 02215 

Abstract 

We show that the argument of Alexander and McTague, that the bcc crys- 
talline structure is favored in those crystallization processes where the first 
order character is not too pronounced, is not correct. We find that any solu- 
tion that satisfies the Alexander-McTague condition is not stable. We inves- 
tigate the implication of this result for nucleation near the pseudo-spinodal 
in near-meanfield systems, 
t Permanent address 

I. INTRODUCTION 

Predicting the symmetry of crystalline phases from first principles and the knowledge 
of the inter-molecular potential has eluded condensed matter physicists for years. One 
promising path was suggested by Alexander and McTague (AM) [jT] who noted that the 
Landau-Ginsburg free energy indicated a particular crystalline form if the first order nature 
of the transition from the liquid was not too pronounced. For systems such as those that 
interact with a Lennard-Jones(LJ) potential the AM argument indicated that the crystalline 
phase would have a triangular or hexagonal structure in d — 2 while the same argument 
suggested that the bcc lattice would be favored in d = 3. This conclusion is independent 
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of the details of the potential as long as, like LJ, it is spherically symmetric and the fluid 
is monatomic. However, simulations of LJ and hard sphere systems as well as experiments 
on metals and rare gasses indicate that these systems freeze into an fee structure. This 
result is also obtained by the application of density functional theories. |||| 

In a related development Klein and Leyvraz(KL) |J used the AM argument to look at 
nucleation of the crystal from the meta-stable liquid in a mean-field system near the spinodal 
and found a hexagonal or triangular structure for the critical droplet in d = 2 and stacked 
triangular or hexagonal planes or a bec structure in d = 3. However, this droplet was 
unstable. In addition, Groh and Mulder(GM) [|7| also looked at a mean-field system near 
a spinodal and, using bifurcation analysis, also noted that the AM argument said nothing 
about the nature of the stable solid but only indicated that a bifurcation to a bec structure 
in d=3 was favored over bifurcations to other periodic structures but the stability of the bec 
structure was beyond the AM or bifurcation analysis. From these results it is still unclear 
what the meaning or accuracy is of the AM argument. It is the purpose of this paper to 
perform some additional analysis of these results. 

In this paper we present an argument indicating that any crystal satisfying the condi- 
tions in AM must be unstable. This argument, which uses the Landau- Ginsburg approach 
employed by AM does not require that the system be mean-field. In addition we use the 
techniques of AM to investigate the nucleation process near the spinodal like singularity in 
systems with long range interactions. We argue that the droplet with the structure indi- 
cated by the AM argument dominates the nucleation process. The fact that this droplet is 
unstable is, unlike the phases proposed by AM, physically reasonable and consistent with 
other results. In addition we investigate the bifurcation analysis of GM and show that the 
additional bifurcations to other symmetries that they find, which are not found by AM are 
not relevant to nucleation. 

The remainder of this paper is structured as follows; In Section II we describe the method 
of AM. In Section III we present the argument for instability. In Sections IV and V we discuss 
the implcations of this argument for nucleation of the crystal from the liquid and in Section 
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VI we summarize and discuss our results. 



II. ALEXANDER-MCTAGUE 

In this section we give a brief description of the argument of AM. |1| We will use a 
slightly different, but equivalent approach. We begin with the Landau- Ginsburg free energy 

F (p) = n I I C{\x - y\)p(x)p{y)dxdy + f f(p(x),T)dx-h f p{x)dx (1) 
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In the above we assume, following AM, that the interaction term is quadratic in p and that 
f(p,T) can be expanded in a power series in p for fixed T. The quantity h is the chemical 
potential. The Euler-Lagrange equation that specifies the equilibrium state is obtained by 
functionally differentiating eq.(0) and setting the derivative equal to zero. 

C(\x - y\)p{y)dy+ ^§p- ~h = (2) 

We assume that there is a constant solution p corresponding to the stable or metastable 
liquid at the values of T and h we are investigating. Defining ip(x) through p(x) = p + ijj(x) 
and inserting this expression in eq.(0) we obtain 

A(\x-y\)^(y)dy + b(h,T)^ 2 (x) = (3) 

Here b(h,T) is the coefficient of the term quadratic in ip(x) arising out of a Taylor series 
expansion of . The dependence of b(h,T) on h and T is through p Q . The term A(\x — 
y\) = C(\x — y\) + b\5(x — y) where b\ is the coefficient of the term linear in ip(x) in the 
Taylor expansion which also depends on h and T through p Q and 5(x — y) is the Dirac delta 
function. Following AM we have truncated the Taylor expansion after the quadratic term 
using the assumption that ip(x) << 1, i.e. the first order nature of the transition is not too 
pronounced. We will see that this can be justified in a self consistent manner. 

Again following AM we consider the system in the neighborhood of a critical point 
associated with the fluid. |TJ The critical point is characterized by A(|fc )|) = 0, where \k Q \ 
is the location of the global minimum of (7(|(A;|) defined in eq.(|l]). Near the critical point 



A(|& |) = e << 1. The hat denotes the Fourier Transform and e = T T J S where T$ is the 
critical temperature. From eq. (|3J) we obtain, after Fourier transform, 

A(\k\)i>(k) + b{h, T) J #(k - k')^(k')dk' = (4) 

There is one solution of eq.(|]) that scales with e i.e. ip(k) ~ e or ip(x) ~ e. The scaling 
justifies neglecting the higher order terms in ip since they will be higher order in e. Such 
scaling solutions, assuming of course that they exist, will have a free energy cost which 
differs from the liquid by terms of order e. What about solutions that do not scale? Unless 
the critical point Ts is on the coexistence curve it seems likely that there would be solutions 
with a lower free energy than the liquid as e — > 0. This existence of other solutions is what 
one would expect if, as we contend, the AM solution is unstable. However, here we are 
simply following the AM argument and do not need to consider the possibility of solutions 
of the Euler-Lagrange equation that do not scale and hence require higher order terms in 
the Landau-Ginsburg free energy then those allowed by AM. The question of solutions for 
nucleation droplets that do not scale with e will be addressed in Section V. 
We now assume that the solution of eq. ([|) is of the form 

^P(x)=J2 c n(h,T)exp(ik n -x) (5) 

n 

— * 

where the k n are reciprocal lattice vectors. If we now take d reciprocal lattice vectors, e.g. 
ki...kd where d is the spatial dimension, to have a magnitude \k a \ then we have a solution 
in which C\...c& » c n for n > d when e << 1; that is for values of h and T near the critical 
point. 

To see this we simply insert eq.(^) into eq.(|3]) to obtain 

oo oo 

c n (h,T)A(\k n \) exp(ik n ■ x) + b(h,T)[J2 c n (h,T) exp(ik n ■ x)f = (6) 

n=l n=l 

We now take Ci(h,T)....Cd(h,T) ~ e and all other c n (h,T) to be higher order in e. We can 
assume A(|A;„|) ~ 1 for n ^ 1, ...,d. It is straightforward to generalize this argument to the 
case of a finite number of degenerate global minima. For e << 1 the modes with magnitude 
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\k \ dominate the solution. In the limit e — ► these modes are the only ones that contribute 
Since these modes dominate for e << 1 they must be such that (Z)n=i ex P(^o,n • %)) 2 has 
the same symmetry as exp(zfc D) „ • x). Here k , n is one of the d reciprocal lattice vectors with 
magnitude \k Q \. Note that the sum contains exponentials and their complex conjugates as 
the sum must be real. Therefore as in AM we must have that the reciprocal lattice vectors 
form equilateral trangles. In two dimensions this generates a triangular or hexagonal lattice 
in real space and in three dimensions either bcc or layered triangular or hexagonal structures. 
This result was also derived via a bifurcation analysis by GM. |7j In the next section we 
investigate the stability of the solutions of eq.(^|). 

III. STABILITY OF THE ALEXANDER-MCTAGUE SOLUTIONS 

To check the stabilty of the solutions of eq.(^) we perform a simple linear stability 
analysis. Writing ip(x) + rj(x), where n(x) is a small arbitrary perturbation and ijj(x) is a 
solution to eq.(|]), and inserting this in eq.(||) we obtain 

A(\x-y\)r}(y)dy + 2b(h,T)-4)(x)ri(x) = $ (rj(x)) (7) 

The function $>(rj(x)) can be expanded in the eigenvectors of the operator on the left hand 
side of eq.(^). If the eigenvalues of the operator are all positive then the solution is linearly 
stable. If there is at least one negative eigenvalue the solution is linearly unstable. We will 
show that there must exist at least one negative eigenvalue. 

The proof makes use of the well known fact that if we have an Hermetian operator © 
defined on a Hilbert space then the average of 0, which we will call Q, which is defined as 

G = <XI&IX> (8) 
< x\x > 

is an upper bound for the lowest eigenvalue of the operator Q. Here \x > is any vector in 
the Hilbert space. The derivation of this result is quite straightforward and can be found 
in Modern Quantum Mechanics by J. J. Sakurai M as well as most elementary books on 
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Quantum Mechanics. Note that the operator in eq.(|7|) is both real and symmetric and hence 
Hermetian. The Hilbert space we use is defined with plane wave normalization. 

Since we are primarily interested in the situation in which the fluid phase with density 
p is stable or metastable we begin with this case first. Since A(k) is the k dependent 
susceptibility of the uniform fluid if the fluid is stable or metastable then A(k) > for all 
\k\. Using if)(x), the solution to eq.@, as the vector \x > in our bound we have that the 
bound B is given by 

= fip*(x)A(\x - y\)jj(y)dxdy + 2b(h, T) J ^(x)^ 2 {x)dx 

J\iP(x)\ 2 dx U 

Since if)(x) is a solution to eq.@ the bound B reduces to 

d_ - SV{x)A(\x-y\)j){y)dxdy 

!\ijj{x)\ 2 dx 1 ' 

Converting to Fourier space by using Parseval's theorem we have 

-B J \4j{k)\ 2 dk = J \^(k)\ 2 A(\k\)dk (11) 

Since is assumed to be positive definite we must have B < 0. Hence, the upper bound 

of the lowest eigenvalue is less than zero which proves the result. Note the fact that ip(x) is 
periodic was never used in the argument which is therefore valid for any solution of eq.([§). 

Although it is not relevant for the AM argument, where the question of the structure of 
the solid phase near a stable or metastable liquid was considered, for the sake of completeness 
we demonstrate that the solutions of eq.(|3]) are unstable when < for some value(s) 

of \k\; that is, when the liquid phase is unstable. Equation |] for the lowest eigenvalue bound 
can be written as 

B [ \w{x)\ 2 dx= [ w*(x)A(\x - y\)w{y)dxdy + 2b(h,T) [ ^(x)w*(x)w(x)dx (12) 



where the operator is given in eq.(|7|). We choose the test funtion w(x) = exp(ik Q ■ x) so 
that 

BV = A{\k \)V + 2b{h,T)C (13) 



where C = J ip(x)dx and V is the system volume. Since k Q is the location of the global 
minimum of A(|A;|) and, by assumption, A(|A;|) < for some \k\, A(\k \) < 0. Returning to 
eq.(|3]) multiplying by b(h,T) and integrating with respect to x yields 

mT)C= ^hixmmi (14) 

V ; A(0) 



Substituting eq.flli]) into eq.(|13D we obtain 

w = A(yy-*P^ (15) 

Vl U A(0) 

Since b(h,T) is real and A(|fc |) is assumed to be negative, if A(0) > then B < 0. If 
A(0) < eq.(0) can no longer be used to generate an equlibrium crystal. This occurs 
because A(x — y) and b(h,T) are functions of p Q which is unstable to spatially constant 
perturbations. 

The one situation left to address is the one in which A(|fc|) is positive for all \k\ except for 
\k \ where it is equal to zero. This is the situation in which the system is at the critical point. 
From eq. (|TT|) we now have B = if ip(k) = S(k — k ). This implies that ip{x) = exp(ik ■ x). 
However, from our discussion in Section II ijj(x) ~ e = A(|fe |) = 0. Therefore the solution 
of eq. (0) in this case describes the infinitely long lived fluctuations at the critical point and 
not a phase. 

^ — # — * 

Note that the condition that we be near a critical point (A(|fc|) = e << 1 for some \k\) 
was necessary to argue that the AM crystals were in fact minima of the free energy. However, 
the critical point condition was not used to show that the solution of eq.(|3]) was unstable. 
In addition, we have not specified the nature of the critical point. When the critical point 
is a spinodal then the AM argument becomes quite useful in determining the nature of 
the critical droplet in the nucleation process that takes place near the pseudo-spinodal in 
near-meanfield systems. We discuss this application in the next section. 

Note also that the same argument that was used to show that the solutions of eq. (§) are 
unstable can be used on any equation of the form 

A(\x-y\)Hy)dy + q(h,T)r(x) = (16) 



where n in any integer greater than or equal to two. 



IV. SPINODAL NUCLEATION OF A CRYSTAL FROM THE LIQUID 

In this section we apply the ideas of AM to the problem of the nucleation of a crystalline 
solid from the liquid near a pseudo-spinodal in a near-meanfield system. First we note that 
spinodals, in the sense that one finds them in the van der Waals theory of liquids or the 
Curie- Weiss theory of magnetic systems, are meanfield(MF) objects. In Ising models this has 
been seen in Monte Carlo simulations |9| and via transfer matrix techniques. [JTTJ A general 



argument using a Ginsburg criterion has also been given by Binder To study MF 



systems in a rigorous way Kac JT2| introduced the idea that, if a system has an interaction 
potential of the form V(|af|) = Ke/d^l) +7 d 0(7|^l) where V re f(\x) is a short range reference 
potential, 7 is a parameter and / •y d (j)('y\x\)dx = D < 00, then in the limit 7 — > the system 
will be MF. This means that for fluids such a potential results in the van der Waals equation 
with the attendant MF critical exponents and spinodals. In magnetic systems the result is 
the Curie- Weiss description. We have assumed that the interaction is spherically symmetric 
for simplicity but that need not be the case. 

In order to generate near-meanfield(NMF) systems we use the approach of Kac however 
we take 7 << 1 but finite. That is we do not take the 7 — > limit. In NMF systems 
there is no true spinodal but the, depending on the interaction range R = 7 _1 the system 
will behave as if there is a spinodal as long as one does not approach the singularity too 
closely. P HTl|1 We will refer to such apparent singularities as pseudo-spinodals. To study 
nucleation in systems with long range interactions undergoing deep quenches we adopt the 
techniques of saddle point evaluation of the partition function [TH] to nucleation near the 



pseudo-spinodal. p]14]| We begin with a Landau- Ginsburg- Wilson Hamiltonian identical to 
F(p) in eq.(|l]) but with one additional requirement. We take the interaction term 

C(\x-y\)= 1 d \{ 1 \x-y\) (17) 

where the A has the properties of the long range Kac potential described above. We will 



take 7 to be small but finite so that we are describing NMF rather than MF systems. Our 
Hamiltonian H is then 

H(p) = \l 1 d k( 1 \x-y\)p(x)p(y)dxdy + £ / b n (h,T)p n (x)dx - h J p(x)dx (18) 

71=1 

where we have made the Taylor series expansion of f(p,T) explicit. The partition function 
in the canonical ensemble is 

Z = J 5pexp(-f3H(p)) (19) 

where (3 = 1/KgT. 

Using R = 7 _1 and assuming that p(r) = p(x/R) the Hamiltonian H in terms of scaled 
lengths becomes 



H(p) = R° 



-J h{\r — ?\)p{r)p{r)drdr + ^5 n (/i,T) J p n {r)dr — h J p{r)dr 



(20) 



n=l 

For R » 1 (i.e. the NMF limit) the partition function in eq.(|T9|) can be evaluated as a 
saddle point integral. The Euler-Lagrange equation that specifies the saddle point will be 
of the same form as eq.(0). Identifying p Q as the constant solution that specifies the density 
of the liquid we write p(r) = p a + ip(r). As in AM we assume that ^(r) is small since we are 
near a pseudo-spinodal critical point. Expanding in tfj(r) we obtain an equation for tfj(r) of 
the same form as eq.@. 

R d [J A(\r - f\)il;(f)df + b(h,T)ip 2 {r)} = (21) 

where A(\f — r*|) = A(|r — r*]) + b\5(r — f") 

To obtain the nucleation droplet we assume a solution of the form 

— > 

^VEaP^'^T) ( 22 ) 

n L 

where k' Q = Rk and L is a length to be determined. We assume at the outset that L » 
|fc | _1 , which will be seen to be true self consistently. As before \k Q \ is the location of 
the global minimum of Since \k \ ~ R, is a slowly varying function of r. 
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Consequently, we can expand ip'(j-) in a gradient expansion about r 1 = r. The first three 
terms in the expansion are 

—f —i ( — * — i\ —*t l — * — v 1 9 — i 

ff iy' I iy> iyi' I tyi' itn ty \ JJ iy" 

= + " V^( Z )|^+ F^^L.)^ (23) 

where the gradient and Laplacian are with respect to j- 

Since A(\f— r*\) — > as |r — f| — > oo, ^ r << 1 for large L and we can truncate the 
series in eq.fp3|) after terms of the second order. Inserting eq. (|22|) and eq. (|23|) into eq.(plf) 
we obtain 

R d [E n exv(ik 0>n -r^A(\k \)^)- 
G(h, T) Enexp(ifc 0jn ■ r)^VV(i) + b (h, T) (En exp(i£ , n ■ r))V 2 (|)] = (24) 

where the Laplacian is with respect to £ and 

G(h,T) = J\r-7\ 2 A{\r-7\)d{f-f*)<0. (25) 

The term involving the gradient is zero as it reduces to a function proportional to 
V^A(|fc|)|^_g and \k \ is the global minimum of Since we are near a spinodal 

A(|A; |) = e << 1. We can assume || that ~ e as long as L ~ e~ 1//2 Hence 



where the scaling of x with respect to both e~ 1//2 and R is explicit. For further discussion 
of this scaling see section V. The argument of AM is now invoked to limit the k ^ n to lie on 
equilateral triangles. The critical or nucleating droplet then has the following form: In the 
interior it is periodic with a triangular or hexagonal structure in two dimensions and a bcc 
or layered hexagonal or triangular plane structure in three dimensions. These structures are 
modulated by an envelope which satisfies the equation 

G(h, T) V 2 ^'(e 1/2 f) + A(\k \)^'(e 1/2 r) + b(h, T)^' 2 (e 1/2 r) = (27) 

where the Laplacian is with respect to r. The boundary conditions are that -?/>'(e 1//2 r*) — ► 
when \ f\ — » oo and that ■^rs'ip'(e 1 ' 2 r) = at |rj = The first condition is simply a statement 
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that the droplet is localized. The second boundary condition assures that the droplet has 
no unphysical kinks at its center. |TJ] Note also that the structure of eq.([27|) implies that 



the solution is a function of \r\. Since the sign of b(h,T) merely sets the sign of ^'(e 1 / 2 ^) 
we lose no generality by assuming that b(h,T) < 0. 

From the discussion in Section III this droplet will be unstable, not just on its surface 
as in classical nucleation [13[], but in its interior. We can see this explicitly by finding the 



eigenvector associated with the negative eigenvalue [jy| of the operator obtained by a linear 
stability analysis about the critical droplet. That is we want the solution of 

J A(\f- 7\)w{r)'d7 - 2\b{h,T)\tP{e 1/2 r)w{r) = Xw{r) (28) 

where ?/;(e 1//2 f) is the solution to eq.fl2~T|) given in eq. (B^) . 
We assume that the eigenvector has the form 

w (f) = exp(ik Q ■ r) W(e l/2 r) (29) 

n 

Employing the same arguments we used to obtain eq.(^) we find that the eigenvector has 
the same interior structure as the critical droplet and an envelope that is the solution of 

G(h,T)V 2 W(e 1/2 r) + eW(e 1/2 r) - 2\b(h,T)\^'(e 1/2 r)W(e 1/2 r) = XW(e 1/2 r) (30) 

where we have set A(|/c D |) = e. 

In d — 1 the solution of eq.(^Tj) can easily be seen by substitution to be 

De 1 / 2 

= h2( m \ ( 31 ) 

cosh {ae 1/z x) 

where a and D are constants. In d = 3 eq. fl2"T|) has been solved numerically |TS| . The solution 
is radially symmetric, has its maximum at the origin and decreases to zero as \r\ — >• oo Hence 
eq.(|3T)D has the form of a Schrodinger equation with a potential V(|f|) given by 

7(M)= k\ko\) -2V5'(6 1/2 H) (32) 

From the discussion above, this potential is a shallow well. The lowest eigenvalue will 
correspond to a bound state and hence will be negative. Moreover the eigenvector corre- 
sponding to the bound state will be centered in the center of the well. From the work of 
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Langer Jl3| we know that the eigenstate that corresponds to the negative eigenvalue is an 
unstable mode of the droplet. The fact that the bound state is centered in the center of the 
well is confirmation that the entire droplet is unstable consistent with the adaptation of the 



AM argument. This form of the eigenvector was also seen numerically in Ising models. |L5 

The dominant contribution to the probability of the occurance of a critical droplet comes 
from inserting the solution to the Euler-Lagrange equation (eq.(U)) into the expression for 
the partition function. The details of this caluclation for liquid-gas and magnetic systems 
can be found in Langers paper for classical nucleation Jl3| and in Unger and Klein for 
spinodal nucleation. The adaptation to liquid-solid nucleation near the spinodal of the 
supercooled liquid of the saddle point part of the calculation is straightforward. || The free 
energy barrier to nucleation then is given by 



F{h,T) = R d J A(\r- ^|)V(e 1/2 f)^(e 1 /V)rffrfr y - \b(h,T)\ J ^(e l/2 r)dr (33) 

where ?/>(e 1//2 r) is the solution to eq.(|2~l"D From eqs. |22| and ^6| we have 

F(h,T) = R d e 3 JA'(\f- F\)\£ n exp{i% 0>n ■ We 1/2 ^)][E„exp(ii% in • r')^{e 1 l 2 r')\drdr' 

- IK^^I/Enexp^-f)] 3 ^ 3 ^)^ (34) 

where A(\f — f*|) = eA'(\f — r'l) and k' on = Rk ^ n 
When \k , n \ ^ 

J exp(z£ 0fn • fWe 1/2 r)dr ~ ^ = e d ' 2 (35) 

where ^ = e -1 / 2 is the correlation length near the spinodal in scaled units. Therefore, the 
dominant contribution to the cubic term in eq.(|34"D will be given by those terms in the sum 
(En e,i °'" T ) 3 which are spatial constants. The same analysis can be done for the quadratic 
term in eq. (|34"D . The dominant contribution will come from the k Q ^ n vectors that appear 
with opposite signs in the sums and result in exponentials of the form exp(/c Di „ ■ (r — r')). 
These terms will be of order £ d and all other terms will be reduced by a factor of as in 
eq-©- 
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The free energy barrier is then 

F{h,T) = MR d e 3 ~ d/2 (36) 
where M is a constant. The probability of a critical droplet is 

P D = g (h, T) exp(-MR d e 3 - d/2 ) (37) 



where g(h,T) is small compared to the exponential. fll"3,14 



The exponential dominates the probability of a nucleation or critical droplet. To get a feel 
for the magnitude of the argument of the exponential we turn to consideration of the so-called 
Ginsburg criterion for the validity of a mean- field treatment. [IT] In MF systems fluctuations 



can be ignored when calculating thermodynamic quantities. |TT| Ginsburg pointed out that 



this implies that the fluctuations of the order parameter must be small compared to its mean 
value. That is 

81 « 1 < 38 > 

where £ is the correlation length as above, xt is the susceptibility and is the order param- 
eter. Using the MF exponents for the spinodal, Xt ~ e" 1 , ~ e and £ ~ -Re -1 / 2 , where we 
have made the R dependence of the correlation length explicit, we have 

« 1 (39) 



R d e 2-d/2 



or 

for MF theory to be a good approximation. If the left hand side of eq. (f40"D is infinite then 
MF theory is exact. Since spinodals are MF objects and pseudo-spinodals affect the physics 
only if the system is NMF we require the condition in eq. ([!(]) to hold if we are to see the 
spinodal nucleation process described above. 

If R » 1 then e can be small and the spinodal can be approached (e << 1) with 
the Ginsburg criterion, R d e 3 ~ d ^ 2 , still valid. For Ising models the spinodal can be seen by 
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measurements of the isothermal susceptibility which will diverge as R — > oo and e — > 0. || 
In the supercooled liquid the static structure factor is known rigorously to diverge at \k\ ^ 0. 
||i~6|j However, in this case the situation is more complicated and the measured divergence is 
suppressed for d > 1. [[H]] 

The implications of these investigations is that the argument of AM, while it does not 
predict the structure of stable crystals, can be adapted to determine the structure of at 
least one kind of nucleation droplet near the pseudo-spinodal in liquids with long range 
interactions. Questions have been raised in the work of GM about the possibility of other 
forms of nucleation droplets near the spinodal. In the next section we will address this 
question. 

V. UNIQUENESS OF CRITICAL DROPLET 

The question of the uniqueness of a critical droplet, given the thermodynamic parameters 
that specify the metastable state, has not been fully resolved. This is true even in phase 
transitions , such as gas-liquid, where there is no spatial symmetry breaking. The resolution 
of this problem is even more difficult when, as in the nucleation of the crystal from the 
liquid, the spatial symmetry changes. To completely answer this question for a given T 
and h, all solutions of the Euler-Lagrange equation resulting from setting the functional 
derivative of H(p) found in eq. (|20|) would have to be found and the free energy cost of each 
solution evaluated. This is a formidable task that has yet to be done either analytically or 
numerically. 

In the simpler liquid-gas nucleation process, although it has not been proven, the idea 
that there is only one saddle point separating the metastable and stable state seems quite 
reasonable. In the liquid-solid transition multiple saddle points with different spatial sym- 
metries seems more plausible. In this regard, Groh and Mulder(GM) 0] have performed a 
bifucation analysis of supercooled liquids near the spinodal. They found that the first order 
bifucation analysis yields only a bcc solution but that the second order analysis results in a 
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solution with an fee symmetry. They also found that within the mean-field Landau theory 
the bee structure has a lower maximum of the free energy than the fee structure. The word 
structure was used rather than phase since, as we have seen in Section III, these structures 
are unstable. Groh and Mulder however, do not consider nucleation. It is the purpose of 
this section to provide an argument, unfortunately not a proof, that as the spinodal is ap- 
proached the bee nucleation droplet found in the previous section dominates the nucleation 
process. 

We begin by noting that it is necessary to have the linear size of the critical or nucleation 
droplet equal to or larger than the correlation length. Particles separated by a scale smaller 
than the correlation length are highly correlated. Therefore, the statistical fluctuations 
that begin the evolution up a saddle point hill should involve regions that are at least the 
correlation length size size. The converse is not true. Critical droplets can be large compared 
to the correlation length since critical droplets are rare compared to the fluctuations that 
set the scale for correlations. It has been the case in all simulations that have looked at 
nucleation that the critical droplet size is either equal to or larger than the correlation length. 
[I8|,r9] This being the case then the critical droplet profile if)(x) can be approximated by the 



form 

ij(x) =5>xp(^ n -r)$(|y|) (41) 

n L 

where the k n are the entire set of reciprocal lattice vectors and L > £. This form of the 
droplet clearly is not exact. We would expect some effect of the fact that the envelope 
will decay to zero as |£| — > oo on the symmetry. However, that will be in the tail of the 
droplet, which is expected to have exponential decay [|1^ ,|6|,|19| (also see eq.(|31~D and will have 
negligible effect on the interior symmetry. Near the spinodal where the correlation length 
£ diverges, the interior of the droplet will be unaffected by the envelope. In addition any 
scaling should also be unaffected by the approximation made at the edge of the droplet. 

With arguments essentially identical to those of Section IV the envelope obeys an equa- 
tion similar to eq . (f24|) . 
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E n exp(ifc n • f)A(|4|)$(£) - \G(h,T)\ E n exp(i£ n • r)V 2 $(£) + 
6(/i,T)(E„exp(z£ n • f)) 2 $ 2 (£) + c(/ i ,T)(E n exp(z4 • f)) 3 $ 3 (|) = (42) 

where the Laplacian is with respect to r and G(h,T), b(h,T) and c(h,T) are independent 
of r and fc. The Laplacian term is of order jp. and L > £ ~ e -1 / 2 . We have restricted our 
considerations to Hamiltonians that include terms up to $ 4 (]r) neglecting higher orders in 
$. These higher order terms can be included with no essential change in the argument. If 
the term 

£exp(i£ n .f)3(|£ n |)->0 (43) 

n 

slower than jn as the spinodal is approached, then the Laplacian term can be ignored. This 
would result, assuming that symmetry constraints could be satisfied, in $(r) being a spatial 
constant. Since nucleation droplets are, almost by definition, localized this result would 
clearly be unsatisfactory. The only way that the term in eq. fl4"3]) can go to zero at all is if the 
sum is limited to terms in which \k n \ = \k \ where k a has the same meaning as in previous 
sections. (See eq.(B^)) 

Therefore, we must have 

J2e*P(ik n ■ r)A(\k n \) = ^exp(ifc ,„ • r)A{\k , n \) ~ e (44) 

n n 

as this is the only possible way that this term can go to zero as e — > 0. Since the terms 
linear in $ must scale the same way this implies that ^ ~ e so that L ~ or L ~ ^. If 
b(h,T) is not zero in eq.(^2|) then $(£) must scale as e. If it were to scale as e x with x > 1 
then the non-linear terms in $ in eq.(|42|) can be ignored relative to the linear terms. This 
results in a solution of the form 

$(I) = ^{ e 1/2 r) = d exp(e 1/2 n ■ r) + C 2 exp(-e 1/2 n • f) (45) 

Lj 

where n is a unit vector whose direction is arbitrary. For the droplet to be localized the 
constant C\ must be set to zero. This results in a violation of the boundary condition that 
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the derivative with respect to \r\ is zero at \r\ = 0. Hence this droplet is physically ruled out. 
If x < 1 then the linear terms in $, including the Laplacian term, can be ignored resulting 
again in a non-local solution. In general nucleation from the liquid to the solid is described 
by Hamilitonians in which b(h, T) ^ and hence the AM argument can be used to fix the 
spatial symmetry in the interior of the critical droplet. Note that b(h,T) is a function of h 
and T and is therefore specified only by the parameters that define the metastable state. 

Suppose that b(h,T) = 0. In this guments similar to the ones above imply that 

$(-jr) ~ e 1 / 2 . The $ 4 term in the Hamiltonian is now relevant and the nucleation barrier 
will now scale as R d e 2 ~ d ^ 2 . In addition a simple extension of AM indicates that the droplet 
symmetry will now be a square lattice in d = 2 and an fee lattice, or stacked layers of square 
lattices, in d = 3. 

It is interesting to ask why, except for a particular symmetry of the Hamiltonian, the 
critical droplet has a bec symmetry in d = 3 when the bifurcation analysis of GM indicates 
that an fee bifurcation is also allowed even when the coefficient (b(h,T)) of the cubic term 
in the Hamiltonian is nonzero. Repeating the GM analysis is beyond the scope of this paper 
however, the key to their construction of the fee bifurcation is the use of lattice vectors for 
which \k\ ^ \k Q \. The argument we used made use of 1) the requirement that the droplet 
be localized, 2)the droplet has no kinks at its center and 3)the assumption that the linear 
droplet size L was greater than or equal to the correlation length £. These conditions taken 
together resulted (see equations |43] and |44|) in a restriction of the reciprocal lattice vectors 
to those for which \k\ = \k Q \. These conditions, while necessary for specifying the nucleation 
problem are not needed in the bifurcation analysis of GM. Hence, even though the fee 
bifurcation is allowed it will not occur in a critical droplet as the spinodal is approached. 

It is important to emphasize that we have restricted our consideration to Hamiltonians 
that are in the same class as that considered by AM. In density functional language this 
implies that all direct correlation functions higher than the pair are zero. Whether the 
inclusion of higher order correlation functions changes the result remains to be seen. 
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VI. SUMMARY AND DISCUSSION 



We have argued that the prediction of a stable bcc phase near a spinodal, or pseudo- 
spinodal, of a supercooled liquid by AM is not correct. The predicted bcc phase is unstable. 
However we have used the AM argument to show that the nucleation or critical droplet near 
a spinodal does have a bcc structure. The same argument that was used to show that the 
AM "crystal" was unstable is also used to show that the critical droplet is unstable. This is 
however, physically reasonable for a critical droplet near a spinodal. 

It should also be noted that droplets like the ones predicted in this paper were found 
in molecular dynamics simulations of supercooled Lennard- Jones fluids. |2(| In addition it 
is known from theoretical studies of the liquid-gas transition that as one moves away from 
the pseudo-spinodal the droplet develops a core that appears almost classical. |~5[ This 
phenomenon is consistent with the results of a density functional calculation of liquid solid 
nucleation by Shen and Oxtoby. |21] 

Many systems of technological importance have long range, or effective long range inter- 



actions. These include polymers ||11|| , neutral plasmas and metals. p2| Since critical droplets 
are at least the size of the correlation length, and all statistical lengths, such as the droplet 
diameter, are measured in units of the interaction range R, classical nucleation, where the 



surface tension is assumed to remain non-zero [O] will be strongly suppressed. In these 
systems, spinodal nucleation will dominate the metastable state decay. Since the structure 
of the interior of a spinodal nucleation droplet is not the stable crystal as would be expected 
in the classical process the evolution of this droplet as it grows will be an important step in 
determining the structure of the stable or metastable crystal. 
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